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We define a general scheme for the evolution of fragmentation functions which re- 
sums soft gluon logarithms in a manner consistent with fixed order evolution. We 
present an explicit example of our approach in which double logarithms are re- 
summed using the Double Logarithmic Approximation. We show that this scheme 
reproduces the Modified Leading Logarithm Approximation in certain limits, and 
find that after using it to fit quark and gluon fragmentation functions to experi- 
mental data, a good description of the data from the largest Xp values to the peak 
region in ^ = ln{l/xp) is obtained. In addition, we develop a treatment of hadron 
mass effects which gives additional improvements at large t;. 



1. Introduction 

The extraction of fragmentation functions (FFs) at large and intermediate 
X from experimental data on single hadron inclusive production has been 
successfully performed ^ using the fixed order (FO) DGLAP ^ evolution 
to next-to-leading order (NLO) However, this formalism fails at small x 
due to unresummed soft gluon logarithms (SGLs), and a different analysis is 
required in which these SGLs are resummed. In this contribution we outline 
a general approach ^ which extends the FO evolution of FFs to smaller 
X values by resumming SGLs. We present explicit results for the case 
where the largest SGLs, being the double logarithms (DLs), are resummed 
using the Double Logarithmic Approximation (DLA) ^'^ in the LO DGLAP 
evolution, and use this scheme to describe, and to fit FFs to, experimental 
data. We also introduce a novel and simple approach to incorporate hadron 
mass effects. 



•Work supported in part by DFG through Grant No. KN 365/3-1 and by BMBF through 
Grant No. 05 HT4GUA/4. 
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2. SGL Resummation in DGLAP Evolution 

The DGLAP equation reads 

^i,(x,<3')=/;i^p(.,<..w>))i,(i,o»), 

where, for brevity, we omit hadron and parton labels. D is a vector con- 
taining the gluon FF Dg and the quark and antiquark FFs Dq and D-g 
respectively, in linear combinations according to the choice of basis, and P 
is the matrix of the splitting functions. We define Os = as/ (27r). We choose 
the simplest basis, consisting of valence quark FFs D~ = Dq ~ Dq, non- 
singlet quark FFs Djvs, and D = {D^,Dg), where for nj quark flavours 
the singlet quark FF is given by _Ds = X^^ii ^q- 

Analytic operations are often simpler in Mellin space, where a function 
f{x) becomes /(w) = /q dxx'^ f(x), since the convolution in x space in Eq. 
(1) becomes the simple product 

-^^D{lo, Q^) = P{uj, as{Q^))D{uj, Q^). (2) 

Consider the formal expansion of P{as) in Ua keeping x (or oj if we want 
to work in Mellin space) fixed, viz. 

oo 

P{as) = J2a^P^''-'^ (3) 

n=l 

where the p("~i) are functions of x (or uj). Equation (3) truncated at some 
chosen (finite) n is known as the FO approach, and, in x space, is not valid 
at small x due to the presence of terms which in the limit a; — > behave like 
(a!^/.T) In^"-"-! x for TO = 1, 2n - 1. Such logarithms arc called SGLs, 
and TO labels their class. As x decreases, these unresummed SGLs will spoil 
the convergence of the FO series for P{x,as) once ln(l/a;) = 0{as ). 
Consequently the evolution of D{x,Q^) will not be valid here, since the 
whole range x < y < I contributes in Eq. (1). 

SGLs are defined to be all those terms of the form a"/a;^"~™ only, 
where m = 1, 2n and labels the class of the SGL, in the expansion about 
w = of Eq. (3) in Mellin space. For to = 1, ...,2n — 1, this definition 
agrees with the form of the SGLs in x space given above, since oj~p = 
-((-l)P/p!) dxx'^il/x)\nP-'^ X for Re(w) > and p > L Such terms 
spoil the convergence of the scries in Eq. (3) as w ^ 0. To construct 
a scheme valid for large and small lo (or x), wc write P in the form P = 
pFO _,_pSGL^ where P^°^ contains only and aU the SGLs in P, so that P^° 
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is completely free of SGLs. Since terms of the type p = {m = 2n), which 
are included in our definition of SGLs, are non-singular, they may therefore 
be left unresummed. Thus we separate pSGL ^^^^ pSGL ^ pSGL _|_ pSGL 
-Pp=o"' which is independent of w, is expanded as a series in Og. On the other 
hand, by summing all SGLs in each class m, Pp^i{w,as) is resummed in 
the form 



oo 



P.f^(-,a.)=E(^) 9m {^), (4) 

m=l 

and truncated for some finite m. The remaining FO contribution 

to P, P^'^{uj,as), is expanded in keeping oj fixed, P^'-*(ti>, Og) = 
Er=i<-P^°'""^H'^)> and truncated for some finite n. Finally, the re- 
sult for P{uj, tts) is inverse Mellin transformed to obtain P{x, Ug), and then 
Eq. (1) is solved exactly. We shall call this the SGL-|-F0-t-F0(5 scheme, 
where "+F0(5" means that the p = terms, which are each proportional to 
6{1 — x) in X space, are left as a FO series in a^. In summary, this scheme 
is the result of resumming all SGLs for which m = 1, 2n — 1 in the form 
of Eq. (4), and treating all remaining terms as in the FO approach. 

In the region for which ttg <C 1 and x is above values for which 
ln(l/x) = 0{as^^^), the SGL+F0(+F0(5) scheme should give a good de- 
scription of the evolution for the following reasons. Firstly, Eq. (4) implies 
Pp>ii^^ <^s) = 7^ J2m=i ("-5 Ina;)" /„ (a^ In^ x). Since Ina; is always 
small, this series is a valid approximation when x is small. On the other 
hand, as a; — > 1 the SGLs for the types p > 1 all vanish, and therefore so 
does each term in the scries. The full contribution from the type p = 
terms is just Pp3j"(.T,a,) = (5(1 — a;) X^^i ^n*^"' ^nd the expansion of 
P^°{x,as) in a„ P^"{x,a,) = J2n=i<P^°^"'^Hx), is finite for ah x. 

3. DLA Improved DGLAP Evolution 

From the DLA ^, if the evolution is rewritten in the form 



dlng^ 



Jx V V V \y 

+ ['^P{y,as{Q'))D(^,Q 



(5) 



then P{x, Qs) is free of DLs. Explicitly, yl = for the DL evolving parts of 
the components D = D~ and D = D^s, while 

/Q 2CV \ 
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for the component D = {D^,Dg). Note that A is a projection operator, 
i.e. it obeys A? = A. P can be determined from P cxphcitly by expanding 
the operator y <ii"Q^ in Eq. (5) in jjj^ and using Eq. (1). 

We will examine what constraint Eq. (5) provides for P. In Mellin space, 
Eq. (5) becomes 

After substituting Eq. (2) into Eq. (7) and dividing out the overall factor 
of D{u), Q'^), we obtain the following constraint on P: 

+ ^dl^) iP-P)+HP-P)P- '^CAasA = 0. (8) 

We now use the fact that P is free of DLs, to obtain an explicit constraint for 
pDL_ £j,g^ make the replacement P = P + P^^ in Eq. (8). Expanding 
Eq. (8) as a series in Us/co keeping as/co^ fixed and extracting the first, 
0{{as/io)^), term gives 

2(pDL)2 ^ ^pDL _ 2CAasA = 0. (9) 

Equation (9) gives two solutions for each component of P. Since P is never 
larger than a 2x2 matrix in the basis consisting of singlet, gluon, non-singlet 
and valence quark FFs, there are four solutions. We choose the solution 

P°'^(a;, as) = ^ (-w + + 16CAa«) , (10) 

since, in the component D = {D^,Dg), the expansion of the result in Eq. 
(10) in as to 0{a^s) keeping w fixed, 

/ n 4Cf „2 I&CfCa \ 

P^\.,as) = r '^^ ""J^ +0{a% (11) 

agrees with the DLs in the literature, while in the components D = D~ 
and D = -Dat^, P^^ — 0. The other possibilities do not give these results 
and/or cannot be expanded in a^, i.e. they are non-perturbative. Equation 
(10) agrees with previous results 
At small CO, Eq. (10) implies that 

Dg,g = (12) 
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obtained after integrating Eq. (2) and neglecting the constant, which is 
valid at large Q. We will use Eq. (12) to partially constrain our choice of 
parameterization at low x in the next section. 

The 0{as) single logarithm (SL) contribution to P is a type p = term 
(see Sec. 2), and is given by 

--'M =(,.:„, _u,:%„ J. m 

Approximating P by agP^^^"' Eq. (7) can be regarded as a generalized 
version of the Modified Leading Logarithm Approximation (MLLA) ^'"^'^ 
equation to include quarks, in the sense that the g component of this latter 

equation for D = (Ds, Dg) when Eq. (12) is invoked is precisely the MLLA 
equation. We therefore conclude Eq. (7) is more complete than the MLLA 
equation. 

P^^iXjUs) is obtained by inverse McUin transform of Eq. (10), which 
yields P^^{x,a,) = Ji (4VCVllni), where Ji{y) is the Bessel 

function of the first kind, given by Ji{y) = ^ Jq dO cos{y sin 9 — 9). 

In the next section we shall consider the DL+L0+L0(5 scheme, in which 
we take P^^l ~ pDL^ ^^^^ ^^^^ pFO (^^j pSCL) ^^^^^ ^^i^^ 

In Fig. 1. we see that Pgg{x, Us) in the DL+LO scheme, which is equal to 
the DL+L0+L0(5 scheme when x ^ 1, interpolates well between its 0{ag) 
approximation in the FO approach at large x and P^J"{x, Ug) at small x (the 
small difference here comes from P^'-'(°)(x) at small x). DL resummation 
clearly makes a large difference to P at small x. 

4. Comparisons with data 

In this section, we study the outcome of fitting FFs evolved in the LO 
DGLAP approach and in our DL+L0+L0(5 scheme to experimental data 
on the normalized differential cross section for light charged hadron pro- 
duction in the process e+e" — > (7, Z) ^ + X, where h is the observed 
hadron and X is anything else. These data exist for different center-of-mass 
(COM) energies and values of Xp = 2p/\/fl. where p is the momentum 
of the observed hadron, which constrain the FFs in the region of x for 
which Xp < X < 1. We fit to data for which ^ < ln(y^/lGeV), where 
^ = ln(l/xp). At LO in the coefficient functions, these data are described 
in terms of the evolved FFs by 
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1 2 3 4 5 6 7 

In (1/x) 



Figure 1. (i) Pgg{x,as) calculated in the DL+LO(+LOi5) scheme, (ii) Pgg{x,as) calcu- 
lated to 0(as) in the FO approach (labelled "LO"), and (iii) P^g^ix, a^) (labelled "DL"). 
as = 0.118/(27r). 



where Qq is the electroweak charge of a quark with flavour q and (Q) is 
the average charge over aU flavours. We wiU take n/ = 5 in afl our calcula- 
tions. Since we sum over hadron charges, we set Dg = Dq . These data de- 
pend on the FFs in the combinations fucix, Qq) — ^ {u{x, Qq) + c{x, Qq)), 
fdsb{x,Ql) = \{d{x,Ql) + s{x,Ql) + h{x,Ql)) and the gluon g{x,Ql). 
For each of these three FFs, wc choose the parameterization 

f{x.,Ql) = Nexp[-cAn^x]x°'{l-x)'^, (15) 

since at intermediate and large x the FF is constrained to behave like 
/(a;, Qo) ~ Nx°'{l — x)^, which is the standard parameterization used in 
global flts at large x, while at small x (where (1 — x)^ « 1) the FF is con- 
strained to behave like lima;^o fix, Qq) = Nexp [— cln^ i — aln i] , which 
for c > is a Gaussian in ln(l/a;), as predicted by the DLA for sufBciently 
large Qq. We use Eq. (12) to remove four free parameters by imposing the 
constraints c„c = Cdsb = Cg and a„c = otdst = ctg. The relation 

Nuc « A^d.b « (16) 

implied by Eq. (12) is not imposed, since we want to describe large x data 
as well. We also fit Aqcd- We choose Q'^ = s, although it is only important 

that the latter two quantities are kept proportional, since the constant of 
proportionality has no effect on the final FF parameters and the description 
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of the data (and therefore the quality of the fit). This implies that there 
will be an overall theoretical error on our fitted values for Aqcd of a factor 
of 0(1). Since all data will be at ^/s > 14 GeV, we choose Qo = 14 GeV. 
The evolution is performed by numerically integrating Eq. (1). 

4.1. Fixed Order Evolution 

We first perform a fit using standard LO DGLAP evolution. We obtain 
Xdf = 3.0, and the results are shown in Fig. 2 and Table 1. The result 
for Aqcd is quite consistent with that of other analyses, at least within 
the theoretical error. It is clear that FO DGLAP evolution fails in the 
description of the peak region. The unphysical negative value of fi for the 
gluon is because the gluon FF is weakly constrained, since it couples to the 
data only through the evolution (see Eq. (14)). 

Table 1. Parameter values for the FFs at Qo = 14 GeV 
pajrameterized as in Eq. (15) from a fit to all data listed 
in the text using DGLAP evolution in the FO approach 
to LO. Aqcd = 388 MeV. 



-...Parameter 


N 









FF 


a 


c 


9 


0.22 


-0.43 


-2.38 


0.25 


u + c 


0.49 


2.30 


[-2.38] 


[0.25] 


d + s + b 


0.37 


1.49 


[-2.38] 


[0.25] 



4.2. Incorporation of Soft Gluon Resummation 

We now redo the previous fit, but now evolving in the DL+LO+LOJ 
scheme. The results are shown in Table 2 and Fig. 3. We obtain = 2.1, 

Tabic 2. Parameter values for the FFs at Qo = 14 
GeV parameterized as in Eq. (15) from a fit to all 
data listed in the text using DGLAP evolution in the 
DL+LO+LO(5 scheme. Aqcd = 801 MeV. 



— -...F'arameter 


N 


P 






FF 


a 


c 


9 


1.60 


5.01 


-2.63 


0.35 


u + c 


0.39 


1.46 


[-2.63] 


[0.35] 


d + s + b 


0.34 


1.49 


[-2.63] 


[0.35] 



Februciry 2, 2008 4:40 Proceedings Trim Size: 9in x 6in 



paper 



20 



15 



• OPAL 202 






A OPAL 172 






• DELPHI 161 






- OPAL 133 






■ OPAL 91 






□ TASSO 44 

■ TASSO 35 

■ TPC 29 




• TASSO U /' 


























'""^ 

i: \ 











Figure 2. Fit to data as described in Table 1. Some of the data sets used for the fit are 

shown, together with their theoretical predictions from the results of the fit. Data to the 
right of the vertical dotted lines have not been used in the fit. Each curve is shifted up 
by 0.8 for clarity. 



15 



• OPAL 202 
A OPAL 172 

4 DELPHI 161 

• OPAL 133 

■ OPAL 91 

□ TASSO 44 
< TASSO 35 

■ TPC 29 

• TASSO 14 




2 3 4 



Figure 3. Fit to data as described in Table 2. 



a significant improvement to the fit above with FO DGLAP evolution. In 
particular, the data arotmd the peak is now much better described. 
We note that had we made the usual DLA (MLLA) choice Q = 
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instead of our choice Q = \/s which is usually employed in analyses using 
the DGLAP equation, we would have obtained half the result Aqcd ~ 800 
MeV. Ng is too large by a factor of about 2 relative to its prediction in Eq. 
(16). However, as noted before, the initial gluon FF is weakly constrained 
in our fits. 



4.3. Incorporation of Hadron Mass Effects 

To attempt to improve the description beyond the peak, we now study 
hadron mass effects, which are important at small Xp. It is helpful to 
work with light cone coordinates, in which any 4- vector V is written in the 

form V = {V+,V-,Yt) with = ^(F" ± V^) and Vt = {V^,V^). 
In the COM frame, the momentum of the electroweak boson takes the 
form q = The momentum of the hadron in the COM frame 

is chosen as = (~^^ "T^v*'^)' ^^^'-^ ^ ~ Ph/l~^ ^® ^^S^* cone 
scaling variable. Therefore the relation between the two scaling variables 
in the presence of hadron mass is Xp = rj ^1 — . As a generalization of 

the massless case, we assume the cross section we have been calculating is 
(da/d^Ms), i.e. ^irj,s) = j]^ ^^(y,s,Q^)D (^Q'), which is related 
to the measured observable (da/dxp){xp, s) via 

do' , ^ 1 da , , . . 

^(xp,s) = — ^^-(ry(x,),s). (17) 

sr;2(xp) 

Although the data is for light charged hadrons, the vast majority of 
particles are pions, so we will assume all particle masses are equal. We 
now perform the DL+L0+L0(5 fit again but with wih included in the list 
of free parameters. We obtain the results in Table 3. The parameters 
are not substantially different to those in Table 2. The result for ruh is 
reasonable for light charged hadrons. We find Xdf = 2.03, i.e. no significant 
improvement to the quality of the fit, and the comparison with data is 
similar to that in Fig. 3. However, treatment of mass effects renders the 
value of Aqcd more reasonable. 



5. Conclusions 

Using the approach in this contribution gives a much better fit to all the 
data than the FO approach and the MLLA ^ do, even if the fit is still not 
in the acceptable range. Further improvement in the large ^ region can 
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Tabic 3. As in Tabic 2, but incorporating mass effects 
in the fit. Aqcd = 399 McV and = 252 McV. 



^~~~-------...Parameter 


N 


/3 






FF 


a 


c 








9 


1.59 


7.80 


-2.65 


0.33 


u + c 


0.62 


1.43 


[-2.65] 


[0.33] 


d + s + b 


0.74 


1.60 


[-2.65] 


[0.33] 



be expected from the inclusion of higher order SGLs. Our scheme ahows 
a determination of quark and gluon FFs over a wider range of data than 
previously achieved, and should be used to extend the NLO global fits of 
FFs to lower Xp values. 
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